Dynamical structure factors of S = 1/2 two-leg spin ladder systems 
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We investigate dynamical properties of S = 1/2 two-leg spin ladder systems. In a strong coupling 
region, an isolated mode appears in the lowest excited states, while in a weak coupling region, an 
isolated mode is reduced and the lowest excited states become a lower bound of the excitation 
continuum. We find in the system with equal intrachain and interchain couplings that due to a 
cyclic four-spin interaction, the distribution of the weights for the dynamical structure factor and 
characteristics of the lowest excited states are strongly influenced. The dynamical properties of two 
systems proposed for SrCu2 03 are also discussed. 

PACS numbers: 75.40.Gb, 75.40.Mg, 78.70.Nx, 75.10.Jm 
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S = 1/2 two- leg spin ladder systems with antiferro- 
magnetic interactions have attracted great attention both 
theoretically and experimentally 0] . Fascinating aspects 
of elementary excitation as well as thermodynamic prop- 
erties have been revealed. Using several theoretical meth- 
ods, it was shown that S = and S = 1 two-triplet 
bound-states exist below the two-triplet continuum in 
addition to the one-triplet excitation |2|, |^, || [| ||, |^, |) . 
Such an S = two-triplet bound-state was identified 
in the optical conductivity measured for an S = 1/2 
two-leg spin ladder material (Ca, La)i4Cu2404i In 
this experiment, the strength of the interchain coupling 
(J_l) and the intrachain coupling (Jy) was estimated as 
J||/J ± - 1 - 1.2 with J± - 1020 - 1100cm- 1 [g [bi|. 

In addition to the interchain and intrachain couplings, 
a cyclic four-spin interaction, which acts among four 
5=1/2 spins forming a plaquette, has been introduced 
to explain experimental findings for cuprate two-leg spin 
ladder systems. The analysis of the one-triplet mode 
observed by inelastic neutron-scattering experiments for 
LagCa8Cu2404i revealed that a cyclic four-spin inter- 
action is necessary to reproduce the observed disper- 
sion relation [O, |l2). For a S = 1/2 two- leg spin lad- 
der material SrCu203, two sets of coupling constants, 
(i) Jj. = 86meV,J|| = 172meV || and (ii) J ± = 
150meV, J|| = 195meV, J cyc = 18meV, Jdiag = 3meV 
Jl4| with J cyc and Jdiag being the coupling constants 
for a cyclic four-spin interaction and a diagonal interac- 
tion, have been proposed to reproduce temperature de- 
pendence of the susceptibility. It seems thus difficult to 
decide the proper model only from temperature depen- 
dence of the susceptibility. Detailed information on dy- 
namical properties and low-lying excitations is desirable 
to discuss characteristics of 5 = 1/2 two-leg spin ladder 
systems. 

In this paper, we calculate the dynamical structure fac- 
tor (DSF) of S = 1/2 two-leg spin ladder systems using 
continued fraction method based on Lanczos algorithm. 
Dynamical properties of S = 1/2 pure two-leg spin ladder 
systems have already been studied in Refs. 15, 3, 16, 17, 
and 4. We perform systematic calculation for DSF and 
characteristics of the lowest excited states, laying stress 
on the effects of a cyclic four-spin interaction. 



Let us first consider the S — 1/2 two- leg spin ladder 
systems described by the following Hamiltonian, 
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where Si t i denotes the S = 1/2 spin operator in the i 
th rung of the 1 — 1,2 chain, and N is the total number 
of the site. The periodic boundary condition is applied 
along the chain. The energy is measured in units of J|| . 
The DSF can be expressed in the form of a continued 
fraction [|l8l as 



S^(q x ,q y ;uj) 
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where l^o) is the eigenfunction corresponding to the 
lowest eigenvalue Eq, q x and q y denote the wave num- 
bers along the leg and the rung, respectively, and z = 
to + ie + E with h = 1. We set e = 3.0 x 10~ 2 . The 
Fourier transform of the spin operator S£ x q is given 

by S£ i9b = (1/VN) exp[i(fej + q v l)]sfj. ^(&, q y ) 
is the static structure factor and C J, (q x ,q y ;uj) is repre- 
sented in the form of the continued fraction. Since rota- 
tional symmetry around the x, y, and z axes remains, we 
only calculate S z (q x ,q y ;uj). 

In Figs. 1(a) and 1(b), we show S z (q x , ir; u>) and 
S z (q x ,0;iv) for N = 28, respectively, in J±/J\\ — 0.5,1, 
and 2. The weight of S z (q x , q y ; uj) is proportional to the 
area of the full circle. The accuracy of the continued 
fraction becomes worse with decreasing J±/J\\ and/or in- 
creasing lo. In Jj_/ Jii = 0.5 the convergences have the rel- 
ative errors about O(10~ 3 ) for lu > 4 and about O(10~ 10 ) 
for lo < 3, whereas in J±/Ju = 2 they have the relative 
errors about O(10~ 10 ) even for u ~ 6. 

In Fig. 1(a), the largest weights in given q x lie in the 
lowest excited states, and the largest weight at q x — -k 
becomes dominant among them with decreasing Jx/Ju. 
The solid line represents the dispersion relation of the 
one-triplet excitation adequate for J±/J\\ ^> 1 0. The 
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are caused by such dressed spinons. 

In Fig. 1(b), the broken line in J±/J\\ = 2, which 
represents the dispersion relation of the S — 1 two-triplet 
bound state obtained for J±/J\\ ^S> 1 @, agrees with the 
numerical results. Also in q y = 0, the largest weights 
in given q x lie in the lowest excited states. As shown in 
Fig. 2(b), the residue in J±/J\\ — 2 decreases drastically 
around q x ~ 0.67T with decreasing q x . This result means 
that the weights spread in many excited states in q x < 
0.67T. Therefore, in Jj_/J\\ = 2, the lowest excited states 
in q x > 0.67T form an isolated mode of the S = 1 two- 
triplet bound state, while those in q x < 0.6-7T become a 
lower bound of the two-triplet continuum. In Jx/Jh < 1, 
the size dependence of the residue becomes noticeable 
except for the region q x ~ 0. Thus, the lowest excited 
states in J±/Jh < 1 are the lower bound of the continuum 
formed by dressed spinons. 

We now investigated the effects of a cyclic four-spin 
interaction described by the following expression, 

N/2 

7~lcyc — Jcyc y~^{4[(<S f i | j ■ Si,j+i)(S2,i+l ■ S2,i) 



FIG. 1: S z (q x , q y ,uj) for S — 1/2 two-leg spin ladder systems 
in (a) q y = n and (b) q y = for TV = 28. Parameters used are 
shown in the figures. The intensity of each pole is proportional 
to the area of the circle. The solid and broken lines represent 
the dispersion relations of the one-triplet excitation B and 
5=1 two-triplet bound-state Q, respectively, adequate for 
J±/J\\ > 1. 



numerical results and the solid line agree with each other 
quantitatively. To discuss characteristics of the lowest 
excited states, we investigate the finite-size effects of 
the poles and their residues of the continued fraction 
C z (q x ,q y ;u>) p9| . As reported in Refs. 19-21, a pole 
which belongs to a continuum tends to have appropriate 
size dependence at least either of its position or of its 
residue. In Fig. 2(a), we show the size dependence of the 
residue for the lowest excited states and that of the pole 
for the lowest and the second lowest excited states. In 
J±/J\\ = 2, no size dependence of C z {q x ,ir;uj) is seen in 
the lowest excited states, while clear size dependence ap- 
pears in the pole for the second lowest excited states. The 
results show that the lowest excited states form an iso- 
lated one-triplet mode below the excitation continuum. 
In Jj_/J|| = 1, the size dependence is slightly seen around 
q x = in the residue. Thus, the lowest excited states in 
q x ~ may be a lower bound of the excitation contin- 
uum, while those in the other wave numbers form an 
isolated mode. In J±/J\\ = 0.5, the residues show the 
size dependence in < q x < 7r, indicating that the low- 
est excited states become a lower bound of the excitation 
continuum. In Jj_/J|| < 1, the elementary excitation may 
be described by two spinons in either chain dressed by a 
weak interchain interaction. The weights in J±/J\\ < 1 
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+(Si,< • S 2 , i+ i) + (S M+1 • S 2 ,i) + -}. (3) 

From inelastic neutron-scattering experiments for 
La6CagCu2404i, the coupling constants were estimated 



as Jj_ = J|| = HOmeV and J, 
we set Jj_/J\\ = 1 and J cyc / J\\ 



eye 



16.5meV |1J. Thus, 
0.1. Combining the 



Hamiltonian (1) with the contribution (3), we calculate 
DSF. The results for N = 28 are shown in Fig. 3. 
In the poles and their residues for q y — n, noticeable 
size dependence comes out in contrast with the case of 
Jx/J\\ = 1 aud Jcyc — shown in Fig. 2. Therefore, in 
q y = 7r the isolated mode is reduced by a cyclic four-spin 
interaction, since it may act as a frustration term: The 
lowest excited states in q x < 0.7tt become a lower bound 
of the excitation continuum, while those in q x > 0.7n 
keep an isolated mode. 

In Fig. 4, the weights of the lowest excited states in 
Jcyc/ J\\ — 0.1 are compared with those in J cyc = 0. In 
q y = 7r, the weights in q x ~ ir are enhanced by J C y C , 
while those in q x < n are suppressed. In q y = 0, on 
the contrary, the weights in q x > 7r/2 are suppressed, 
while those in q x < ir/2 are enhanced. It is apparent 
that the distribution of the weights of DSF is drastically 
rearranged by a small cyclic four-spin interaction. 

We next investigate dynamical properties of the sys- 
tems described by the parameter set (ii) for SrCu203 
mentioned previously [H. We add the contribution of a 
diagonal interaction described as Jdiag XX^i.i • S2A+1 + 
S\ ; i + i ■ S2s) to the expressions (1) and (3). The results 
for N = 28 are shown in Fig. 5. Note that the re- 
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FIG. 2: The size dependence of the poles and their residues corresponding to the results in Fig. 1. 
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FIG. 3: S z (q x ,q y ,ui), and the size dependence of the poles 
and their residues in J±/J\\ = 1 and J cyc / J\\ = 0.1 for N = 28. 



suits for the parameter set (i) {Jx/Jh = 0.5) [[13) have 
been already presented in Figs. 1 and 2. In q y = it and 
0, noticeable size dependence appears both in the poles 
and the residues. Therefore, the lowest excited states 
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FIG. 4: The weights of the lowest excited states for J cyc / Jy 
0.1 and in J±/J« = 1. 



in < q x < 7r are the lower bounds of the excitation 
continuum as in the case of the parameter set (i). 

In Fig. 6, the dispersion relations and the weights of 
the lowest excited states for both parameter sets (i) and 
(ii) are compared. In q y = n, the dispersion relation 
for the parameter sets (ii) becomes natter in q x < 0.7it, 
whereas the weights for both parameter sets are almost 
the same. The results are qualitatively the same as those 
for N — 24 p2]| . In q y = 0, the dispersion relation for 
the parameter set (ii) becomes flatter with suppressed 
weights in q x > ir/2. The difference in the distribution of 
the weights between two parameter sets may be difficult 
to observe experimentally. On the other hand, it has been 
shown that the effects of a fiat dispersion relation of the 
lowest excited states appear in phonon-assisted optical 
absorption spectrums [|, |[ [hJ. To decide the proper 
model for SrCu203, such experiments may be effective. 
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In summary, we have investigated dynamical proper- 
ties of S = 1/2 two- leg spin ladder systems. It has been 
shown that the largest weights in given q x lie in the lowest 
excited states. The distribution of the weights is influ- 
enced by a cyclic four-spin interaction. We have further 
shown that due to a cyclic four-spin interaction, the low- 
est excited states for q y = tt in LagCasCu2404i become 
a lower bound of the excitation continuum in q x < 0.77T, 
while those in q x > 0.7ir may form an isolated mode. We 
have also discussed the difference in dynamical properties 
between two systems proposed to describe SrCu203. 



FIG. 5: S z (q x ,q y ,ui), and the size dependence of the poles 
and their residues in J±/J\\ = 0.770, J cyc / J\\ = 0.0923, and 
Jdiag/J\\ = 0.0154 for N = 28. Those values correspond to 
the parameter set (ii) for SrCu2C>3 H]. 
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FIG. 6: The dispersion relations and the weights of the lowest 
excited states in the parameter sets (i) and (ii) for SrCu2 03. 



N.H. would like to thank K. Okunishi for useful com- 
ments. Our computational programs are based on TIT- 
PACK Ver. 2 by H. Nishimori. Numerical computation 
was partly carried out at the Yukawa Institute Computer 
Facility, Kyoto University, and the Supercomputer Cen- 
ter, the Institute for Solid State Physics, University of 
Tokyo. This work was partly supported by a Grant-in- 
Aid for Scientific Research from the Ministry of Educa- 
tion, Culture, Sports, Science, and Technology, Japan. 



[1] E. Dagotto, Rep. Prog. Phys. 62, 1525 (1999), and ref- 
erences therein. 

[2] O. P. Sushkov and V. N. Kotov, Phys. Rev. Lett. 81, 
1941 (1998). 

[3] K. Damle and S. Sachdev, Phys. Rev. B 57, 8307 (1998). 
[4] V. N. Kotov, O. P. Sushkov and R. Eder, Phys. Rev. B 



59, 6266 (1999). 
[5] C. Jurecka and W. Brenig, Phys. Rev. B 61, 14307 
(2000). 

[6] S. Trebst, H. Monien, C. J. Hamer, Z. Weihong, and R. 

R. P. Singh, Phys. Rev. Lett. 85, 4373 (2000). 
[7] W. Zheng, C. J. Hamer, R. R. P. Singh, S. Trebst, and 



5 



H. Monien, Phys. Rev. B 63, 144410 (2001). 
[8] C. Knetter, K. P. Schmidt, M. Griininger, and G. S. [15] 

Uhrig, Phys. Rev. Lett. 87, 167204 (2001). 
[9] M. Windt, M. Griieninger, T. Nunner, C. Knetter, K. P. [16] 
Schmidt, G. S. Uhrig, T. Kopp, A. Freimuth, U. Ammer- 
ahl, B. Biiechner, and A. Revcolevschi, Phys. Rev. Lett. [17] 
87, 127002 (2001). [18] 

[10] M. Griieninger, M. Windt, T. Nunner, C. Knetter, K. P. 

Schmidt, G. Uhrig, T. Kopp, A. Fre imuth. U. Ammerah l. [19] 
B. Biiechner, and A. Revcolevschi, :ond-mat /0109524 . [20] 

[11] M. Matsuda, K. Katsumata, R. S. Eccleston, S. Brehmer, 

and H.-J. Mikeska, Phys. Rev. B 62, 8903 (2000). [21] 

[12] S. Brehmer, H.-J. Mikeska, and M. Miiller, N. Nagaosa, [22] 
and S. Uchida , Phys. Rev. B 60, 329 (1999). 

[13] D. C. Johnson, Phys. Rev. B 54, 13009 (1996). 

[14] Y. Mizuno, T. Tohyama, and S. Maekawa, J. Low. Temp. 



Phys.117, 389 (1999). 

T. Barnes, E. Dagotto, J. Riera, and E. S. Swanson, 
Phys. Rev. B 47, 3196 (1993). 

D. B. Yang and W. C. Haxton, Phys. Rev. B 57, 10603 
(1998). 

R. Eder, Phys. Rev. B 57, 12832 (1998). 

E. R. Gagliano and C. A. Balseiro, Phys. Rev. Lett. 59, 
2999 (1987). 

M. Takahashi, Phys. Rev. B 50, 3045 (1994). 

H. Yokoyama and Y. Saiga, J. Phys. Soc. Jpn 66, 3617 

(1997). 

N. Haga and S. Suga, Phys. Rev. B 65, 014414 (2002). 
Y. Mizuno, T. Tohyama, and S. Maekawa, J. Phys. 
Chem. Solids 62, 273 (2001). 



